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We investigate a Friedmann universe filled with a tachyon scalar field, which behaves as dustlike 
matter in the past, while it is able to accelerate the expansion rate of the universe at late times. 
The comparison with type la supernovae (SNIa) data allows for evolutions driving the universe into 
• a Big Brake. Some of the evolutions leading to a Big Brake exhibit a large variation of the equation 

' of state parameter at low redshifts which is potentially observable with future data though hardly 

detectable with present SNIa data. The soft Big Brake singularity occurs at finite values of the 
scale factor, vanishing energy density and Hubble parameter, but diverging deceleration and infinite 
pressure. We show that the geodesies can be continued through the Big Brake and that our model 
universe will recollapse eventually in a Big Crunch. Although the time to the Big Brake strongly 
depends on the present values of the tachyonic field and of its time derivative, the time from the Big 
Brake to the Big Crunch represents a kind of invariant timescale for all field parameters allowed by 
SNIa. 
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INTRODUCTION 



' Dark energy (DE) models aim to explain the accelerated expansion rate of the universe at late times. This phe- 
nomenon was originally discovered using supernovae la (SNIa) data [l[ and has since then been confirmed by many 
other observations (see e.g. (H and references therein). Still, the nature of dark energy, and the precise physical 
mechanism producing the accelerated expansion remains to date an outstanding mystery for cosmologists and for 
theoretical physicists. 

The ACDM model, based on a cosmological constant and cold dark matter, appears to be in good agreement with 
most of the present observational data on large cosmological scales. However, this model has well known theoretical 
problems Q and it also encounters difficulties in explaining some of the data on the scales of structures and even on 
very large scales, like peculiar flows (see e.g. [J]). 

Alternatives to the ACDM model are dark energy models with a time varying equation of state (EoS) parameter 
w and these are not yet excluded by data. In this context, many scalar field models have been considered, either 
5—1 ' minimally coupled with a standard kinetic term, or more complicated ones, e.g. Dirac-Born-Infeld (DBI) models with 
kinetic terms involving a square root 5] . Interest in DBI type models was revived in the framework of string theory, 
where the respective scalar fields arc called tachyons [6], 7]. These models are possible dark energy candidates, as they 
can be interpreted as perfect fluids with a sufficiently negative pressure in order to produce the late-time accelerated 
expansion. 

There is a large arbitrariness in the choice of the potential for tachyonic cosmological models. In Ref |8| (to be 
referred henceforth as I) a specific tachyon potential, containing trigonometric functions, was considered. This model 
turns out to be surprisingly rich, in that it admits a large variety of cosmological evolutions depending on the choice of 
initial conditions. Thus, in I two interesting properties were found. First, for positive values of the model parameter 
k, the sign of the pressure can change during evolution. Second, while under certain initial conditions the universe will 
expand indefinitely towards a de Sitter attractor, under different initial conditions, after a long period of accelerated 
expansion the pressure becomes positive and the acceleration turns into deceleration. Accordingly, the tachyon field 
will drive the universe towards a new type of cosmological singularity, the Big Brake, characterized by a sudden stop of 
the cosmic expansion. At this singularity, the universe has a finite size, a vanishing Hubble parameter and an infinite 
negative acceleration. In contrast to this dynamical picture, similar evolutions were analyzed earlier from a purely 
kinematical standpoint and named sudden future singularities Q . As already stressed earlier for some tachyon models 
0, the addition of dust is a non-trivial problem as is also the behavior of the model at the level of perturbations. So 
our model cannot be viewed yet as a fully viable cosmological scenario but rather as a toy model that could lead to 
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a viable one after suitable improvements. In a recent paper [10(1 (to be referred henceforth as II) we have confronted 
our tachyon cosmological model with SNIa data [ll[ (see also [l2j|). The strategy was the following: for fixed values of 
the model parameter k, we scanned the pairs of present values of the tachyon field and of its time derivative (points in 
phase space) and we propagated them backwards in time, comparing the corresponding luminosity distance - redshift 
curves with the observational data from SNIa. Then, those pairs of values which appeared to be compatible with the 
data, were chosen as initial conditions for the future cosmological evolution. Though the constraints imposed by the 
data were severe, both evolutions took place: one very similar to ACDM and ending in an exponential (de Sitter) 
expansion; another with a tachyonic crossing where the pressure turns positive from negative, ending in a Big Brake. 
It was found that a larger value of the model parameter k enhances the probability to evolve into a Big Brake. For 
a set of initial conditions favored by the SNIa data, we have also computed in II the time to the tachyonic crossing, 
and the Big Brake, respectively These time-scales were found to be comparable with the present age of the Universe. 

The purpose of the present paper is twofold. First we propose to shed more light on the evolution of the tachyon 
field in the distant and in the more recent past; and second, to explore in detail what happens when the universe 
reaches the Big Brake. As this singularity is a soft one, with only the second derivative of the scale factor diverging, 
it is expected that it may be possible for geodesic observers to cross the singularity. Indeed, the traversability of a 
rather generic class of sudden future singularities by causal geodesies was put in evidence in [l3|. Strings can also 
pass through [HI ]. 

In Section II we consider the late-time evolution of the tachyon field, its energy density p, pressure p and EoS 
parameter (barotropic index) w defined as p = wp. In particular, we investigate whether some observable signature 
today may point towards a Big Brake singularity in the future. 

In Section III we discuss the Big Brake singularity, both in terms of curvature characteristics and by analyzing the 
geodesic deviation equation. In Section IV we discuss what happens to the tachyon universe after the Big Brake. 
Finally, we summarize our results with some comments in the Concluding Remarks. 



II. TACHYON SCALAR FIELD COSMOLOGY 



First, we briefly give the basic equations of tachyon cosmology. The Lagrangian of a tachyon field is 



L = -V(T) y /l-g^T li T l/ , (1) 

where V(T) is some suitable tachyon potential. A homogeneous tachyon field T(t) in a Friedmann-Lemaitre- Robertson- 
Walker (FLRW) universe with metric 

ds 2 = -clt 2 + a 2 {t) [dr 2 + r 2 (d9 2 + sin 2 9 dip 2 )] , (2) 

can be thought of as an ideal (isotropic) comoving perfect fluid with energy density p given by 

V(T) . . 

P = i , 3 

Vl-T 2 

and pressure p given by 



p = —V{T)yl - f 2 , (4) 
where a dot denotes the derivative with respect to cosmic time t. The Friedmann equation is then 



rl 8nG V(T) 



H 2 = -- . y '. , (5) 
3 Vl-T 2 



while the equation of motion for the tachyon field T reads 



l-s 2 V 

where 



s Vt 

—j + 3Hs + = , (6) 



s = T. 



(7) 
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Here we consider the following potential 



V(T) 



1- (1 + fc) cos 2 



:VW+k)T 



(8) 



where fc and A are free model parameters. 

^From the present values To and so of the phase space variables T and of its time derivative s — T we found 
convenient in II to introduce the parameters tjq and Xq (denoted Wq in II) defined as 



ya = cos 



; VA (1 + fc) To 



(9) 



1 



X 



1 



The Hubble parameter H as a function of the redshift z is expressed as 



H\z)=H£ n T ,o 



Po 



with Ot o 



p po - , which can be computed in principle as follows 



P(£) 



exp 



, l + w(z') 



3 / cfe' 

o 1 + 



(10) 



(11) 



(12) 



where w = ^ can be obtained from ([3]) and (|4]). 

As mentioned in the Introduction, we consider here a model containing only the tachyon field T. Hence, with respect 
to the expansion rate, the EoS parameter w of T should be compared to what is usually called w e g = wde^de, for a 
universe filled with a dark energy component and dust-like matter. In particular, for ACDM we have w c g = — Oa> 

Finally it is instructive to write down the second Friedmann equation 



a 
a 



AttG 



(13) 



The Big Brake corresponds to vanishing energy density p and infinite (positive) pressure p. Hence d = from ([5]) 
while a = — oo from (|13[) at the Big Brake (whence the name). It is reached for finite time and finite value of the 
scale factor. 



A. Evolution of the system 

We consider the evolution of the trajectories of the model compatible with the supernovae data [ll| at the 1-a 
level. To this purpose we first display on Fig [T] the behavior of the distance modulus and of the luminosity distance 
as functions of redshift z, for 4 different models, all of them fitting within 1-a accuracy the data (actually, the curves 
have the best fits in their respective model classes). We recall that the distance modulus fi is defined as 

^ = 51Ogl0 ^c 1+25 ' (W) 
where di J (z) is the luminosity distance which, for a flat Friedmann universe, is given by (c = 1) 

^> = < 1 + 2 >f < 15 > 

For these specific evolutions we also show the normalized dimensionless energy density p , pressure —2— = w p 

and EoS parameter w (being w e s for the ACDM model) as function of the redshift (both in the past and in the future) 
on Figs [13] 
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FIG. 1: (Color online) The distance modulus (left panel) and dimensionless version (with c = 1) of the luminosity distance 
(right panel) in recent cosmological times. The models shown here are ACDM (black); the tachyonic model with k — —0.4 
evolving into de Sitter (dS, blue); tachyonic models with k = 0.4 evolving into de Sitter (green) and into a Big Brake (BB, 
red) . The curves are all in very good agreement with the SNIa data, as they pass close to the local minimum of the (respective 
regions of the) 1-a domains selected by supernovae (see Figs 2 and 3 of II.) 




FIG. 2: (Color online) The evolution in the recent past of the universe and in the future of the normalized energy density p 
(left panel) and of the normalized pressure p (right panel) is shown for four models. One of the models leads to a Big Brake 
singularity in the future while the other three models shown tend asymptotically to a de Sitter space. 



Four curves appear on eac h g raph on Figs [TJ21 The black curves refer to the ACDM model (with value of 12a, o 
taken from WMAP analysis |l5[)- The dark matter component evolves here as (1 + z) 3 . The other three curves are 
for the toy model containing only the tachyon field. The three curves differ in the model parameter k, and in the 
initial data xq, yo- All three curves pass close to the local minimum of the respective 1-a domains selected by type 
la supernovae. They are characterized by the model parameter k = —0.4 (blue) and k = 0.4, respectively. For the 
latter we have picked up both types of allowed evolutions, one going into de Sitter (green), the other ending in a Big 
Brake (red). Note also that nowadays we have u> = w(z = 0) € [—0.8, —0.6] in all the evolutions displayed. These 
values are similar to ACDM where w e ff,o = — ^a,o ~ —0.74. As expected from viable cosmological evolutions, the 
parameter w approaches zero in the past, corresponding to dust-like behavior at early times. For all parameters, the 
pressure remains very slightly negative for the whole positive z range plotted, thus w stays below zero (although it 
gets very close to it). Under the assumption that this behavior is unchanged when dust is added, we expect to have a 
model where dark energy (here the tachyon field) remains nonnegligible in the early stages of the universe, with 52<r 
remaining roughly constant (actually slightly increasing) during the matter era before it would start dominating at 
late times. A similar behavior can also be achieved in some scalar-tensor DE models, see e.g. [lfjj . 
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k=0.4, y =-0.8, x =0.725, x 2 =195.77, BB k=0.4, y =-0.8, x =0.725, x 2 =195.77, BB 

k=0.4, y =-0.05, x =0.860, £=\ 95.77, dS k=0.4, y =-0.05, x =0.860, x 2 =1 95.77, dS 

k=-0.4, y =0.3, x =0.815, x 2 =195.39, dS k=-0.4, y =0.3, x =0.815, x 2 =195.39, dS 

SJ A =0.726, £! m =0.274, x 2 =1 95.64, dS n A =0.726, Q m =0.274, * 2 «1 95.64, dS 



FIG. 3: (Color online) The evolution of the EoS parameter w (w e s for the ACDM model) in the recent past and in the future 
is shown for the same four models, at late times from z = 2.5 to z = — 1 (left panel); and in the distant past from z = 100 on 
(right panel). It is seen that the tachyon field behaves essentially as dustlike matter at high redshifts. 



B. Observable signature of the Big Brake? 



Though the Big Brake is a singularity taking place in the future, it is interesting to investigate whether models 
leading to a Big Brake exhibit some signature in our present-day universe. As one can see from Fig SJ this is indeed 
the case for some of the universes with a Big Brake in the future which have a characteristic behavior of the EoS 
parameter w in our past: a "dip" at low redshifts. This happens when the Big Brake is not too far in the future, 
in other words when the final redshift is substantially larger than —1. It is then interesting to investigate whether 
such a behavior can be detected. Actually large variations of w at low redshifts can easily hide in luminosity-distance 
curves. Therefore, we expect that it is very difficult to observe this peculiar behavior with present SNIa data. Let us 
consider why this is the case in more details. 

In our model we have only one component, the tachyon field, so what we call here w is actually what should be called 
w e ff when more fluids are present. However, it would be easy to extend our analysis to the case when a dust-like 
component is also present. The EoS parameter w has a minimum at some small redshift z m i n with z m i n ~ 0.2 in the 
most favorable cases. In principle this can be tested and it is straightforward to derive the following equality 



(1 "T" Zmin) 



din/ 



dz 



(16) 



as well as the inequalities 



(l + z)—L < -i 

dz 

. .din/ 
dz 



z < z„ 



z > z n 



(17) 
(18) 



where the quantity / is easily expressed in terms of the luminosity-distance, viz. 



In 



1 



(19) 



a prime standing for the derivative with respect to z. 

Hence it is seen that the (in)equalities (|16 p .p7 p .(|18 p imply a condition on the third derivative of the luminosity 
distance dj}. It is clear that even small uncertainties on di,(z) can lead to large uncertainties on its derivatives 
thereby rendering the (in)equalities ([T6l) .(fT7 |) . (|18p ineffective. In addition, as we have to apply these conditions on low 
redshifts, these involve third order contributions in z to di,(z) which are inevitably small. Obviously, observational 
uncertainties are presently too high in order to make use of (|TB)) . p7|) . (|18l) with existing SNIa data. 
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FIG. 4: (Color online) The evolution of the EoS parameter w in the recent past and in the future for model parameters leading 
to a Big Brake and in the l-cr domain of supernova data. All evolutions have a dip (when, as z decreases, the decrease in w 
turns into an increase), some of them already in the recent past. 

Still, we have tried to see whether a standard x 2 analysis making use only of SNIa data in the range < z < z m in 
could differentiate models with and without dip. This means that we assume a priori that the model with the dip at 
Zmin is the correct one and that we investigate whether a simple statistical analysis of this relevant part of the data 
could hint at its presence. As could be expected, even in this case we find no statistical evidence for the detection 
of a dip with the present SNIa data. Note that even the Constitution dataset [l?} contains only 147 SN data at 
redshifts z < 0.2. Though more refined statistical tools should clearly be used (see e.g. [Hj]), it is quite obvious that 
the present data do not allow for an unambiguous detection. Note that in some models a characteristic smoother 
variation of w can take place on a larger range of redshifts (see e.g. [2(|) which should be easier to detect. Models 
involving a very large variation of w at extremely low redshifts z <C 1 were considered in [19j and it was found 
that they could escape all high precision measurements. In our case, variations, though not as large, are located at 
higher redshifts. We conjecture that future SNIa surveys, like e.g. the Large Synoptic Survey Telescope (LSST) and 
Wide Field InfraRcd Survey Telescope (WFIRST) containing many more supernovae and reducing significantly the 
systematic and statistical errors could allow for such a detection. Future surveys like Euclid involving weak-lensing 
are also promising in this respect. We believe this could be an interesting scientific goal for these surveys especially if 
peculiar models with a large variation of their EoS parameter at low, but not too low, redshifts, like some of our Big 
Brake models with a dip, are in good agreement with observations and are theoretically motivated candidates. 

III. THE BIG BRAKE SINGULARITY 
A. Curvature 

The 3-spaces with t =const have vanishing Riemann curvature 

{3) n abcd = o . (20) 

The 4-dimensional Riemann curvature tensor has therefore but few nonvanishing independent components: 

Rtrtr = —da , 

Rttptip = Rtete sin 2 9 — —aar 2 sin 2 9 , 
Rnprip = R r 0r9 sin 2 9 = a 2 a 2 r 2 sin 2 9 , 

Re v e v = a 2 a 2 r 4 sin 2 9 , (21) 

and the corresponding components arising from symmetry. Remarkably, all components which diverge at the Big 
Brake are of the type Rtata- Therefore, the singularity arises in the mixed spatio-temporal components. 
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B. Geodesic deviation 



The geodesic deviation equation along the integral curves of u — d/dt (which are geodesies with afhne parameter 
t) is 

u a = -R a cbdV b u c u d , (22) 

where r] b is the deviation vector separating neighboring geodesies, chosen to satisfy rfUf, = 0. In the coordinate system 
([2]) we obtain 

u a = -B? fU Ti b <xa, (23) 

which at the Big Brake diverges as — oo. Therefore, when approaching the Big Brake, the tidal forces manifest 
themselves as an infinite braking force stopping the further increase of the separation of geodesies. This can be also 
seen from the behavior of the velocity 

v a = u b V b rf oc H , (24) 

which at the Big Brake vanishes. Immediately after, the negative acceleration will cause the geodesies to approach 
each other. Therefore a contraction phase will follow: everything that has reached the Big Brake will bounce back. 

We conclude the section with the remark that despite the singularity of the geometry (the second derivative a of the 
scale factor diverges at t = ies), its soft character (d stays regular) assures that a continuation of the evolution is still 
possible in the following sense. We indeed need to know a to follow the evolution of the spacetime but we only need 
to know a to follow the evolution of free particles. This means that despite not being able to continue the evolution 
of the geometry in a direct way, we can univocally continue the individual world lines of freely falling test particles 
(geodesies), each of these being perfectly regular at t = tgg. The singularity is not experienced by any individual 
freely falling particle, but makes itself felt only through the equation of geodesic deviation, which at t — tsB indicates 
that the expansion of the geodesic congruence turns negative from positive. 

Once the particles have gone through the Big Brake, we can again start to evolve the geometry itself, thus following 
the further evolution of the universe beyond the singularity. As will be shown in the following section, the tachyonic 
universe will evolve along similar trajectories to those starting from a Big Bang, but in the opposite direction (with 
s — > — s), arriving therefore into a Big Crunch. 



IV. FROM THE BIG BRAKE TO THE BIG CRUNCH 
A. How the universe crosses the Big Brake singularity 



To understand how the crossing of the Big Brake singularity takes place, and what is going on after the crossing, 
it is convenient to refer to the phase portrait of the model at some positive value of the parameter k. This portrait 
was drawn in Paper I and we reproduce it here, see Fig. 5. The accessible phase space of the model consists of a 
rectangle (T3 < T < T4, —1 < s < 1) and four stripes. The values 

2 1 

13 = — , =arccos , . 

3^/(1 + fe)A VTTk 

T 4 = ; 2 [ 7T - arccos ) (25) 

3y/(l + k)A V VTTk) 

are those for which the potential ((5) vanishes. Inside the rectangle the dynamics of the system is described by the 
Lagrangian (JlJ with the potential ((HJ , while in the strips the Lagrangian is given by 



L = W{T) y /g^T ll T v -\, 



(26) 



with the potential 



W{T) = 



A 



sin 



x/AO + fc) T 



'(1 + k) cos 2 



VA(l + fe)T 



1 . 



(27) 
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FIG. 5: (Color online) Phase portrait evolution for k > (k = 0.44) 



Consider a trajectory entering the left lower strip through the point Q' having coordinates (T = T3, s = — 1). The 
analysis of the equations of motion, carried out in Paper I, has shown that the universe encounters a Big Brake (BB) 
singularity after a finite time. This singularity is characterized by some value of the tachyon field T BB , of the time 
t BB and of the value of the cosmological radius a BB . These values are found numerically up to normalization, as was 
done in II. 

The equation of motion ([6]) for the expanding universe in the left lower strip can be written as 

1/4 



3s(s 2 - 1)VA I (k + 1) cos 2 = 



sm 



3^A(l+fc)T 



3y/A(l + k)(s 2 - 1) 3y/A(l + k)T 



(k + 1) cos : 



3^A(l+fc)T 



+ (fc-l) 



(28) 



1 



^From the analysis of this equation we find that, approaching the Big Brake singularity in the lower left strip of the 
phase diagram, the tachyon field T, its time derivative s, the cosmological radius a, its time derivative a and the 
Hubble variable H behave respectively as 



T = T> 



BB 



\3W(T BB ] 



1/3 



(tBB-t) 1 / 3 , 



(29) 



( 4 

\81W(T BB ] 



1/3 



(t B B - t) 



_,\-2/3 



(30) 



a — O-BB — ~^ a BB 1 



f9W 2 {T BB ) 



1/3 



(tBB - tf /3 , 



(31) 



a = a BB 



/ 9W 2 (T BB ) 



1/3 



{tBB - t) 



- ^1/3 



(32) 



H 



9W 2 (T BB ) 



1/3 



(t B B - t) 



- ^V3 



(33) 



9 



To arrive to formulae (l29l) (|33j) we have used the following strategy. Assume that in the neighborhood of the Big 
Brake singularity the tachyon field behaves as 

T = T BB + A(t BB -t) a , (34) 
where A and a are some real parameters to be determined. Then, s behaves as 

s = -aA{t BB -t) a - 1 , (35) 

while its time derivative is 

s = a{u-l)A(t BB -t) a - 2 . (36) 

A simple calculation shows that the first "friction" term, proportional to the Hubble variable in the right-hand side 
of Eq. ([28]) , has the behavior 

s 5/2 ~(tss-i) 5(Q ~ 1)/2 , (37) 

which is stronger than the corresponding behavior of the second potential term in the right-hand side of Eq. 
which is 



,s 2 



(tBB-t) 2 ^. (38) 



This means that the term s in the left-hind side of Eq. ([28]) should have the same asymptotic as the friction term in 
the right-hand side of the same equation and, hence 



which gives immediately 



( 39 ) 



a = i. (40) 



Comparing the coefficients of the leading terms in Eq. ([28]) we find that 



A -(fmsi)- (41) 

Thus, we arrive at Eq. ([29)) . Eq. (|30l) follows right away. Using the Friedmann equation we obtain the value of the 
Hubble parameter (Eq. (|3"3")0 . which, in turn, gives formulae (13"TT) and (|32p for the cosmological radius and for its time 
derivative. 

The expressions (p29"|) - (|3"3")) can be continued in the region where t > t BB , which amounts to crossing the Big Brake 
singularity. Only the expression for s is singular at t = t BB , but this singularity is integrable and not dangerous. 

Upon reaching the Big Brake, it is impossible for the system to stop there because the infinite deceleration eventually 
leads to the decrease of the scale factor. This is because after the Big Brake crossing the time derivative of the 
cosmological radius ([3"2"j) and of the Hubble variable ([3"3"|) change their signs. The expansion is then followed by a 
contraction. 

Corresponding to given initial conditions, we can find numerically the values of T BB , t BB and a BB (see Paper 
II). Then, in order to see what happens after the Big Brake crossing, we can choose as initial conditions for the 
"after-Big-Brake-contraction phase" some value t = t BB + e and the corresponding expressions for T,s,H,a and a 
following from relations (|29l) - (|33[) . and integrate numerically the equations of motion, thus arriving eventually to a 
Big Crunch singularity. 



B. What is going on after the Big Brake crossing ? 



After the Big Brake crossing the universe has a negative value of the variable s, less than —1. This means that 
its evolution should end in a finite period of time. Remember that the universe is now squeezing. The equation of 
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motion for s then looks as follows: 



3s(s 2 - 1) VA / (k + 1) cos 2 3 ^ A( 2 1+fc)T - 1 



1/4 



. 3jA(l+k)T 



3 v /A(l + fc)(s 2 - 1) 3JA{l + k)T 

— cot — 

2 2 



(fc + i)cos 2 VApiz: + ( fc _i) 

( fc + 1)c0S 2 VAW _! 



(42) 



In principle, the evolution of the universe can either end at the vertical line T = at some value of s or at the 
horizontal line s = — 1 at some value of T. One can find the corresponding points on the phase diagram by direct 
analysis of the system of equations of motion. However, such an analysis is rather cumbersome. Thus it is convenient 
to use some results of the analysis of the trajectories for the expanding universe given in Paper I. Begin by writing 
down the equation for the trajectories describing the expanding universe in the phase space (T, s), eliminating the 
time parameter t: 

3(l-s 2 )VA / l-(fc + l)cos 2 ^p y T \ 1/4 

dT sk ^p y i- s 2 



3 y/XpTk) i - s 2 mt / 3^HTTk-)T \ (k + 1) cos 2 3 V^ T + (fc_ i) 

2 3 C ° V 2 / 1 - (* + 1) cos 2 

This equation is valid in both the rectangle of the phase diagram and in the strips. In I we have considered it in the 
upper left strip and saw there that the trajectories of the phase diagram can have their beginning only in the points 



(T = 0, s = yj *±1), (T = 0, s = VkTT) or [T = T*,s = 1), < T* < T 3 . Now, it is easy to see from Eq. @ that 

the simultaneous change of sign of the Hubble parameter a/ a and of the time derivative of the tachyon field leaves 
this equation invariant. This means that the trajectories describing the expansion from the Big Bang singularity 
in the upper left strip are symmetrical reflections with respect to the axis s — of the trajectories describing the 
contraction towards the Big Crunch singularity in the lower left strip. Thus, equation (I43|) written above describes 
also the trajectories of the contracting universe in the lower left strip. In turn, this implies that all the contracting 



trajectories can only end at the points (T = 0, s = (T = 0,s= -VF+1) or (T = T*, s = -1), < T* < T 3 . 

We are now in a position to analyze the behavior of all these trajectories, describing the contracting universe, 
using the results obtained for the corresponding trajectories, born in the left upper strip, describing the expanding 
universe. First, all the contracting trajectories encountering the Big Crunch singularity at the point (T = T*, s = — 1), 
where < < T3 and the unique trajectory ending in the point (T = 0, s = —yk + 1) enter the lower left strip 
from the rectangle of the phase diagram through the corner (T = T^^s = —1) without arriving from the Big Brake 
singularity (indeed, they originate from the repelling de Sitter node). These trajectories are the time- reversed of the 
corresponding expanding trajectories, having their origin at the points (T = T*, s = +1) and of the unique trajectory 
originating at the point (T = 0, s = +\fk + 1), which enter into the rectangle through the point P (see Fig. 5) and 
which do not undergo any change of the expansion regime. 

Thus, the only point where the trajectories coming from the crossing of the Big Brake singularity can end to is 

(T = 0,s = —yj ^-^-)- Now recall (see Paper I) that the trajectories born at (T = 0, s = \J -^r") behave at the 
beginning of their evolution as 



k + 1 2(i-fci 

±Z_ + DT-W, (44) 

where the parameter D can take any real value. Among these, those with a sufficiently large positive value of D, 
say, D sep < D grow without limit and do not achieve a maximal value of the variable s. Instead, they approach 
asymptotically to the vertical line T — Tbb, s —> +00, thus encountering a Big Brake shortly after the Big Bang (such 
a possibility was overlooked in I). Instead, the trajectories, for which D < D sep achieve some maximal value of the 
variable s after which they turn down and enter the rectangle of the phase diagram through P. Such trajectories were 
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TABLE I: Key moments in the evolution of the tachyon universes for k — 0.2. Columns (1) and (2) report different pairs 
of values of the magnitudes yo and xo which are compatible with the supernovae data within la confidence level. Columns 
(3), (4), (5) and (6) report respectively the times t», t BB and tsc elapsing from the present to the tachyonic crossing, to the 
attainment of the Big Brake and to the later attainment of the Big Crunch, and the time lapse between the Big Brake and the 
subsequent Big Crunch. (The values of t», t BB and tsc have been calculated assuming for the Hubble parameter the value 
H = 73 km/s/Mpc.) 



yo 


Xo 


U (W 9 yrs) 


t BB (lOVs) 


tsc (W 9 yrs) 


{tBC-t BB ) {W 9 yrs) 


-0.90 


0.635 


0.334 


1.042 


1.412 


0.198 


-0.85 


0.845 


2.377 


3.093 


3.300 


0.207 


-0.85 


0.860 


2.438 


3.146 


3.352 


0.206 


-0.85 


0.875 


2.505 


3.206 


3.410 


0.204 


-0.80 


0.890 


6.237 


6.927 


7.135 


0.206 


-0.80 


0.905 


6.663 


7.348 


7.554 


0.206 


-0.80 


0.920 


7.197 


7.877 


8.082 


0.205 



described in detail in Paper I. The trajectory characterized by the critical value of the parameter D — D sep plays the 
role of separatrix between these two sets of the evolutions. 

The trajectories approaching the Big Crunch at the point (T = 0, s = — y^p) behave as 



1 „ 2(i-*o 
DT^+^~ 



(45) 

Those with D > D sep are the evolutions which underwent the Big Brake crossing in the left lower strip. 

It is interesting to study the properties of the special cosmological evolution mentioned above (D = D sep ) which 
separates in the upper left strip the subset of trajectories attaining a maximum value of s and then entering the 
rectangle at point P from the subset of those trajectories for which s is not bounded above and which encounter the 
Big Brake already in the left upper strip. This separatrix is composed by two branches, one in the upper left strip 
and a symmetrical one in the lower left strip both having the vertical line T = T3 as asymptote. It encounters the 
Big Brake singularity at T3 at some time moment t = t BB . Now, analyzing Eq. (|28[) for this trajectory, we find that 
in the left neighborhood of t — t BB the dynamical variables behave as follows: 



T = T 3 - A a (t BB -t) 



- *W 



(46) 



s = -A a (t BB 



t) 



-5/7 



(47) 



H = 



U7y/Ak(k+ 1) 

4Ao 



1/4 



BB 



(48) 



where 



Aa= (§) 1/6 A- 5/12 (l + fc)" 5/12 fc 1/4 



(49) 



The analysis leading to these formulae is analogous to the one which led us to formulae (|!?9"|) - (pi5)) . 

If we continue (|4l) |) - (|4"9"|) beyond t = t BB we see that the expansion turns into a contraction, the tachyon field T 
starting to decrease, while its time derivative jumps at t = t BB from an infinite positive value to an infinite negative 
one. Thus, &tt>t BB the universe finds itself in the lower left strip. It leaves the Big Brake along the asymptotic line 

T = T 3 , s = —00, and eventually attains the Big Crunch singularity at the point (T = 0, s = — J ^-)- In the lower 

left strip the separatrix separates the trajectories exiting from the Big Brake singularity at some value T BB < T3 from 
those which enter the strip from the rectangle, attain a minimum value of s and end in the Big Crunch. 

In summary, the expansion phase of the universe, which originated from a Big Bang, stops at t = t BB and turns 
there into a phase of contraction leading eventually to a Big Crunch. 
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TABLE II: As in Table U for k = 0.4. 



yu 




t* (l0 9 iyrs) 


( BR ( 10 9 vrs) 

L LJ LJ V if I 


tnr (10 9 vrs) 


(t Rr - tun) (lO 9 vrs) 


—0.80 


0.710 


0.836 


1.644 


1.933 


0.289 


—0.80 


0.725 


0.841 


1.629 


1.915 


0.286 


—0.80 


0.740 


0.847 


1.616 


1.900 


0.284 


—0.75 


0.815 


2.153 


2.952 


3.247 


0.295 


—0.75 


0.830 


2.195 


2.983 


3.277 


0.294 


—0.75 


0.845 


2.242 


3.020 


3.312 


0.292 


—0.70 


0.845 


3.845 


4.635 


4.932 


0.297 


—0.70 


0.860 


3.964 


4.746 


5.043 


0.297 


-0.70 


0.875 


4.097 


4.871 


5.168 


0.297 


—0.70 


0.890 


4.247 


5.015 


5.310 


0.295 


-0.65 


0.860 


6.182 


6.959 


7.259 


0.300 


—0.65 


0.875 


6.473 


7.243 


7.540 


0.297 


-0.65 


0.890 


6.808 


7.573 


7.870 


0.297 


-0.65 


0.905 


7.204 


7.963 


8.259 


0.296 


-0.60 


0.875 


10.253 


11.016 


11.314 


0.298 


-0.60 


0.890 


11.108 


11.866 


12.163 


0.297 


-0.60 


0.905 


12.203 


12.956 


13.251 


0.295 


-0.55 


0.875 


19.517 


20.274 


20.570 


0.296 


-0.55 


0.890 


25.030 


25.782 


26.077 


0.295 



C. The time lapse from the Big Brake to the Big Crunch 

We give here the Tables I JlIII which, for a selection of values of the model parameter k, report the times elapsing 
from today to the tachyonic crossing, to the attainment of the Big Brake and to the later attainment of the Big 
Crunch for some trajectories which are compatible with SNIa data within a 1-tr confidence level. 

As we can see, although the time to the tachyonic crossing and to the Big Brake both depend strongly on the initial 
data within the l-cr domain selected by supernovae, the time interval between the Big Brake and the Big Crunch do 
not show the same dependence. Nevertheless, this interval exhibits a slight increase with the model parameter k. 

V. CONCLUSIONS 

As shown in II there are tachyon cosmologies described by Eqs.(IT|)- p^)) which are compatible with the supernovae 
data and which are subject to a Big Brake in the future. In this paper we have addressed some questions about this 
model: how it behaves in the distant past, whether these model universes can produce observational signatures today 
and whether they can be continued beyond the Big Brake singularity. 

Having in mind the eventual construction of fully viable models, we are comforted by the fact that the tachyonic 
field has a (quasi) dust-like behavior in the past regardless of its future evolution in the parameter range allowed by 
supernovae. Present supernovae data can hardly discriminate between the evolutions going into a de Sitter phase and 
those leading to a Big Brake. However, we emphasize that the large variation of the EoS parameter at low redshifts 
occurring in some of the evolutions leading to a Big Brake might be detectable with future high precision data. 

The main result of the paper is the study of the cosmological evolutions going into a Big Brake. These evolutions 
can be extended beyond the Big Brake, despite the geometry becoming singular, which apparently forbids such a 
continuation. However, this singularity is a soft one, as only the second derivative a of the scale factor diverges at 
t = tsB, while a does not. We need to know a to follow the evolution of the spacetime but we only need to know 
a to follow the evolution of free particles. This means that we can not continue the evolution of the geometry, but 
we can univocally continue the individual world lines of freely falling test particles (geodesies), each of these being 
perfectly regular at t = tss- The singularity is not experienced by any individual freely falling particle, but makes 
itself felt only through the equation of geodesic deviation, which at t = tsB indicates that the expansion of the 
geodesic congruence turns negative from positive. 
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TABLE III: As in Table [TJ for k = 0.6. 



yu 




U (Vfivrs) 


( BR C 10 9 vrs) 

L LJ LJ V if I 


tnr (10 9 vrs) 

V if I 


(tor - tnn) (lO 9 vrs) 


—0.75 


0.665 


0.548 


1.369 


1.693 


0.324 


—0.70 


0.755 


1.434 


2.289 


2.624 


0.335 


—0.70 


0.770 


1.451 


2.289 


2.623 


0.334 


—0.70 


0.785 


1.469 


2.292 


2.625 


0.333 


—0.70 


0.800 


1.489 


2.299 


2.628 


0.329 


—0.65 


0.815 


2.561 


3.401 


3.740 


0.339 


—0.65 


0.830 


2.614 


3.443 


3.782 


0.339 


—0.65 


0.845 


2.671 


3.490 


3.827 


0.337 


-0.60 


0.830 


3.854 


4.692 


5.036 


0.344 


—0.60 


0.845 


3.960 


4.788 


5.131 


0.343 


-0.60 


0.860 


4.077 


4.897 


5.237 


0.340 


—0.60 


0.875 


4.206 


5.018 


5.359 


0.341 


—0.55 


0.845 


5.510 


6.336 


6.682 


0.346 


-0.55 


0.860 


5.711 


6.530 


6.875 


0.345 


—0.55 


0.875 


5.937 


6.749 


7.091 


0.342 


—0.55 


0.890 


6.194 


6.999 


7.339 


0.340 


—0.50 


0.845 


7.460 


8.281 


8.629 


0.348 


—0.50 


0.860 


7.803 


8.617 


8.963 


0.346 


-0.50 


0.875 


8.193 


9.002 


9.345 


0.343 


—0.50 


0.890 


8.645 


9.447 


9.790 


0.343 


—0.45 


0.860 


10.668 


11.478 


11.823 


0.345 


-0.45 


0.875 


11.370 


12.174 


12.519 


0.345 


—0.45 


0.890 


12.208 


13.008 


13.349 


0.341 


-0.45 


0.905 


13.237 


14.033 


14.373 


0.340 


-0.40 


0.860 


15.044 


15.851 


16.196 


0.345 


-0.40 


0.875 


16.471 


17.273 


17.615 


0.342 


-0.40 


0.890 


18.313 


19.110 


19.453 


0.343 


-0.40 


0.905 


20.838 


21.631 


21.972 


0.341 


-0.35 


0.860 


23.487 


24.291 


24.635 


0.344 


-0.35 


0.875 


27.874 


28.674 


29.016 


0.342 


-0.35 


0.890 


36.194 


36.989 


37.328 


0.339 


-0.30 


0.845 


43.469 


44.276 


44.621 


0.345 



Since the geometry can be reconstructed by the knowledge of each of its geodesies, the evolution of the universe 
does not stop at the Big Brake. Once the particles have gone through the Big Brake, they will determine the geometry 
anew and we can start to evolve the universe beyond the singularity. A phase of contraction follows, leading eventually 
to a Big Crunch. 

We have analytically and numerically analyzed the evolution of the tachyonic universe from the Big Brake to the 
Big Crunch. Quite remarkably, the numerical study showed that although the time to the tachyonic crossing and 
to the Big Brake both depend strongly on the initial data chosen from the I- a domain selected by supernovae, the 
time intervals between the Big Brake and the Big Crunch do not exhibit the same dependence and they only slightly 
depend on the model parameter k (within a factor of two) . This seems to provide an invariant timescale for the class 
of tachyonic scalar cosmologies considered, presumably related to the fact that some information (the behavior of the 
higher derivatives of the scale factor) is lost while passing through the Big Brake. 
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